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Abstract 

The Parisi formula for the free energy in the Sherrington-Kirkpatrick and mixed 
p-spin models for even p > 2 was proved in the seminal work of Michel Talagrand, 
|24j . In this paper we prove the Parisi formula for general mixed p-spin models which 
also include p-spin interactions for odd p. Most of the ideas used in the paper are well 
known and can now be combined following a recent proof of the Parisi ultrametricity 
conjecture in [17J. 
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1 Introduction and main result. 

The formula for the free energy in the Sherrington-Kirkpatrick model [21j was famously 
discovered by G. Parisi in [18\, \i9\ using the approach that combined a replica trick with 
a very special choice of the replica matrix. It was later understood in [9], [10] that the 
special form of the replica matrix conjectured by Parisi corresponded to a number of physical 
properties of the Gibbs measure of the model, one of them being the ultrametricity of its 
support. The Parisi formula for the free energy in the Sherrington-Kirkpatrick and mixed 
p-spin models was proved by M. Talagrand in |24j| following the discovery of the replica 
symmetry breaking interpolation scheme by F. Guerra in [8]. However, for technical reasons 
only the case of p-spin interactions for even p > 2 was considered. Using the main result in 
[T7] which yields that under a small perturbation of the Hamiltonian the support of the Gibbs 
measure in these models is indeed asymptotically ultrametric, we prove the Parisi formula 
for general mixed p-spin models that include odd p-spin interactions as well. Let > 1. Let 
us consider Gaussian processes H]si^p{cr) for p > 1 indexed by cr e Sat = { — 1, +1}^, called 
pure p-spin Hamiltonians, 

^nM) = ^ip-i)/2 9h,...,i,(^ii ■ ■ ■ (^i,: (1-1) 

i<n,...,jp<Af 
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where random variables (s'ii,...,^^) are standard Gaussian independent for all p > 1 and all 
{ii, . . . ,ip). Let us define a mixed p-spin Hamiltonian as their linear combination 

Hn{ct) = Y,(^pHnA^) (1-2) 
p>i 

with coefficients {/3p) that decrease fast enough, for example, Ylp>i '^^l^p < This technical 
condition is sufficient to ensure that the process is well defined when the sum includes 
infinitely many terms. The covariance of the Gaussian process Hj^^a) is easy to compute and 
is given by a function of the normalized scalar product, called overlap, i?i,2 = N''^ '^i<N ^I'^l 
of spin configurations cr^ and <t^, 

^HN{<T^)HN{a^) = m{Ri,2), (1.3) 

where ^{x) — X]p>i/^p^^- Given A; > 1, let us consider two sequences of parameters, 

= mo < rrii < . . . < ruk-i < ruk < 1 (1.4) 

and 

= go < ?i < ■ ■ ■ < gfc < Qk+i = 1, (1-5) 

which will be denoted by m and q and consider independent Gaussian random variables 
{zj)o<j<k with variances = ^'(g^+i) — i'{qj)- We define 

Xfe+i=logch Zj and = — log E/ exp m/X^+i (1.6) 

0<j<k ' 

recursively for / < k, where E/ denotes the expectation in the r.v. {zj)j>i. When mi — this 
means Xi = E;X;+i. Let us denote 9{q) = q^'{q) — i{q) and define 

Vk{m,q)=Xo{m,q)-^ J2 m,{e{qj+,) - e{q,)) . (L7) 

i<j<fe 

Then the following theorem holds. 
Theorem 1. (The Parisi formula) We have, 

lim ^Elog expHN{cr) = MVk{m,q), (1.8) 
where the infimum is taken over all k, m and q as above. 

The quantity in the limit on the left hand side is called the free energy of the model and the 
infimum on the right hand side is the famous Parisi formula. One can include the external 
field term in the model but for simplicity of notation we will omit it. 
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2 Proof. 



Most of the ideas of the proof are well known and available in different places in the literature. 
Under various formulations of the ultrametricity conjecture, one can find arguments that 
contain many of the same ideas in [3] and [I5l in the case of models with only even Jo- 
spin interactions, and a sketch of the proof of the general case in Section 15.3 in [25]. The 
ingredient that was missing is the main result in [17] which also allows us to handle the case 
of the general mixed p-spin models. 

The Ghirlanda-Guerra identities. A central role in the proof is played by the Ghirlanda- 
Guerra identities [7J that are utilized in two distinct ways. First, they yield positivity of 
the overlap via Talagrand's positivity principle which allows us to obtain the upper bound 
using Guerra's replica symmetry breaking interpolation scheme and, second, they imply 
ultrametricity of the overlap array using the main result in [17] which allows us to identify 
the asymptotic Gibbs measures that appear in the proof of the lower bound based on the 
Aizenman- Sims- Starr scheme [1]. Let us consider a perturbation Hamiltonian 

i/rV) = Ar-V8^2-%ff^,^(a), (2.1) 

where if^p(cr) are independent copies of the p-spin Hamiltonians in (11. ip and (xp)p>i are 
i.i.d. random variables uniform on an interval of length one, for example, [1,2]. Replacing 
Hn with H]\i + H^'^^ in (II. 8p obviously does not affect the limit since the perturbation term 
is of a smaller order. However, adding this perturbation term regularizes the Gibbs measure 
in the following way. Let be the Gibbs measure on Sjv corresponding to the Hamiltonian 

Hn + K^^\ 

eMHN{<r) + Hl-\cT)) 
Gn[o-) = , {2.2} 

where Z]^ = X^o-eSiv ^^P(-^^(''') + ^N^^^i^)) denote by (■) the average with respect to 
the product Gibbs measure Gfj^. Let (cr')i>i be an i.i.d. sequence of replicas sampled from 
On and denote by 

^M' = ^E^'^'' (2.3) 

i<N 

the normalized scalar product, or overlap, of cr' and cr''. Given p > 1, n > 2 and a bounded 
measurable function / of the overlaps (-R/,z')/,z'<n on n replicas, let 

^{f,n,p) = \E,{fRl„,^,) - ^E,{f)E,{Ry -lf2^MKi)[ (2-4) 

1=2 

where denotes the expectation with respect to all Gaussian random variables for a fixed 
uniform sequence {xp)p>i. Then, the Ghirlanda-Guerra identities can be stated as follows. 

Proposition 1. For any p > l,n > 2 and a bounded function f of the overlaps {Ri^v)i^v<n, 

lim E,0(/,n,p) = 0, (2.5) 

where E^- is the expectation with respect to {xp)p>i. 
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The proof of this result is well known and we refer to Chapter 12 in |25j for details. We will not 
be using these identities directly for measure Gn but for other Gibbs measures with a slightly 
modified Hamiltonian Hn{(t) since it is well known that the proof of the identities is robust 
to such modifications and depends mostly on the form of the perturbation Hamiltonian f l2.ip . 
It is interesting to note that once we finish the proof of Theorem [H the argument in [14J will 
immediately imply that fl2.5p holds in a strong sense without the perturbation Hamiltonian 
for all p > 1 such that /3p 7^ in (OD- 

Guerra's replica symmetry breaking bound. In the case when p-spin interactions for 
odd p > 3 are not present in (11.21) . the inequality < in (11. 8p was proved by F. Guerra in [8J by 
inventing the replica symmetry breaking interpolation scheme. The fact that this inequality 
holds even in the presence of odd p-spin interactions was observed by M. Talagrand in |23] 
and we will only briefly remind the main idea, which is to write down Guerra's interpolation 
scheme in terms of the Ruelle probability cascades [20] (Poisson-Dirichlet cascades in the 
terminology of [25]) and force the overlap to be positive along the interpolation by adding 
the perturbation term (12. ID . Given k > 1, the Ruelle probability cascades are defined as 

(i) a random probability measure (wo)^^^^ 011 via some explicit construction involving 
Poisson-Dirichlet point processes PD{() with intensity measures x~^~''dx for ( G (0, 1) and 

(ii) a Gaussian process {za)aeN>' with covariance 'Kzaiz^i = ^'(q'aiAaO where 

A a'^ = min{/ > 1 : aj ^ af} if 7^ and A a'^ = k + 1 ii = a'^ 
(see Chapter 14 in [25] for details). For < t < 1 we define an interpolating Hamiltonian 

HN,t{fy, a) = VtHN^cr) + y/l-t ^ Za^i^i (2.6) 

i<N 

where {za,i)aem are independent copies of {za)aem for ? > 1 and let 

ipit) = j^ElogY,w^exp{H^,t{a,a) + H'f;'\(T)). (2.7) 

a,cr 

If we define the Gibbs measure on S^v x N'^ by 

Fi{(cr,a)} ~ Waexp(HN^tio-,a) + H^''\cr)) 
then a straightforward calculation using Gaussian integration by parts gives 

ip'it) = 9(1) + ^E{9iq^.^^2))^^ - ^E(e(i?i,2) - Ri,2C{qa^Aa-) + e{q^i^^2))^^ (2.8) 

where {■)rt is the Gibbs average with respect to Ff^. When ^(x) = ^p>i /S^a;^ does not 
contain terms for odd p > 3, ^ is convex on [—1, 1], which implies that the last term in (12. 8p 
is negative and dropping this term and integrating the corresponding inequality for < t < 1 
we obtain an upper bound on the free energy in (II. 8p . The fact that the representation of 
this upper bound in terms of the Ruelle probability cascades coincides with the formula in 
(II. 7p is well known and is explained in great detail in Chapter 14 in [25]. If the terms for 
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odd p > 3 are present, the function ^ is only convex on [0, 1] but the argument still works 
if we know that -R1.2 is non-negative with high probability under EFf^. This is where the 
perturbation term in (12. 7p comes into play to ensure that the Ghirlanda-Guerra identities 
hold along the interpolation and, as a consequence, to ensure the positivity of the overlap 
via Talagrand's positivity principle (see Section 12.3 in [25]). In fact, an observation in [TT] 
shows that the perturbation term Hf^^* forces the positivity of the overlap uniformly over 
all measures on S^r in the following sense. If given a measure i^iv on Sjy we define a random 
probability measure z>jv on S^v by the change of density dv^i^cr) ~ exp Hf^^\(T) (iz/jv(cr) then 
Theorem 1 in [IT] implies that for any e > 0, 

lim supEz>^2(/2^2 < s) = 0. (2.9) 

Using this for the marginal on T.^ of the Gibbs measure 7t{(cr, a)} ~ Wq, exp H]\f^t{(^, a) 
on Sjv X implies that the remainder term in (12. 8p is asymptotically non-negative and we 
can proceed as in the case of even p-spin interactions. 

The Aizenman-Sims-Starr scheme. The proof of the lower bound is done in several steps 
but it begins with the Aizenman-Sims-Starr scheme p]. Let us consider the Hamiltonian 
^jvl^r) = Ep>i /3p^^,p(o-) where 

^nA"^) = (N + l)(p-m 9H,...,p(^h (2.10) 

^ ' l<ii,...,ip<N 

Let Gj^ and (■)_ denote the Gibbs measure and its average corresponding to the Hamiltonian 
+ Hf^^^ and let z{cr) and y{fT) be two Gaussian processes on S^v with covariances 

Ez{a')z{a') = ^(i^i.a), Ey{(T')y{a') = 6{R,,2) (2.11) 

independent of each other and all other random variables. Then the Aizenman-Sims-Starr 
scheme in [1] yields the following (see e.g. Section 15.8 in [25]). 

Proposition 2. The limit in U.8\) is bounded from below by 

liminf(Elog(chz(<T)) - Elog(expy(cr)) ). (2.12) 

The only difference here is that we included the perturbation term Hf^^^((T) but, since it 
is of a smaller order, one can easily check that it does not affect the computation leading 
to this representation. Below, we will express the limit (12.12^ in terms of some asymptotic 
Gibbs measure that satisfies the exact form of the Ghirlanda-Guerra identities but, in order 
to do so, we first need to show that Propositions [D and [2] also hold with non-random choices 
of the sequence x = {xp)p>i (depending on A^) rather than on average over x. We mentioned 
above that the proof of the Ghirlanda-Guerra identities is robust to modifications of the 
Hamiltonian and, in particular, they hold for the Gibbs measure Gjj so that if in ( 12. 4p 
we replace (■) by (■)_ then (12. 5p still holds. Let us consider a collection 

= i (/,'^,P) : P > 1,^ > 2, / is a monomial of {Ri,v)i^v<n \- 
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Since this a countable family we can enumerate it, {{fj,nj,pj))j>i, and define a function 

<Pt = M^) = J2^~'<P{fj,nj,Pj) (2.13) 

which depends on the variables in a; = {xp)p>i. Since each monomial |/| < 1, we can see 
from the definition (12. 4p that |0(/, n,p)| < 2 and, therefore, the Ghirlanda-Guerra identities 
fl^ imply that 0. Let 

A = A(x) = Eg\og{ch z{cr)) - E,, log(exp y(cr)) (2.14) 

where, again, denotes the expectation with respect to all Gaussian random variables for 
a fixed x. We will need the following simple lemma. 

Lemma 1. We can find x = {xp)p>i such that 

Mx) < c(E,0^)^/2 and A(x) < E,A + c(E,0^)i/^ (2.15) 
where c is a constant that depends only on the function ^. 

Proof. If we denote by E^ and Ey the expectations with respect to {z{(t)) and {y{a)) then 
(12. lip and Jensen's inequality imply 

< E<;log(chz(<T)) < Eglog(E^chz((T)) = ^'(l)/2 

and 

< E3log(expy(o-))_ < Eg log(E,expy(o-))_ = ^^(l)/2 

and, therefore, the function — c < A(x) < c. Given e > 0, let = {x : A(x) < Ea;A + e}. 
Then, E^A > (E^A + e)Px(fi') - cV^iVt) and ^^iVt) > e/(E^A + e + c) > ejc. On the other 
hand, by Chebyshev's inequality, ^xij^T < ^) > 1 — ^x4>tI^ and if e/c > E^0jr/e then 
^ l~l < 7^ 0- Taking e = c(Ea;0jr)^/^ implies that we can find a point x that satisfies 
both inequalities in (I2.15p . 

□ 

For each A^, let us choose = (x^)p>i that satisfies (I2.15P and, since E^0jr — )■ 0, we get 
lim0j-(x^) = O and liminf E^.A > liminf A(x^). (2.16) 

Let us redefine the Hamiltonian H^^^ and the Gibbs measure by fixing parameters 
X = x^ and, since the measure now depends only on the Gaussian randomness, we will 
write E instead of E^. By fl2.16p . Proposition |2] still holds for this redefined measure G]^ and, 
recalling ( I2.13p . 

1 1 " 

ft 11/ 

1=2 

for all p > 1, n > 2 and all monomials / of {Ri,i')i,i'<n- 
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Asymptotic Gibbs' measures. Next, we will define an asymptotic analogue of the Gibbs 
measure and represent the limit (12.121) in terms of this measure. Let (cr')^>i be an i.i.d. sample 
from Gjj and let = {Rj^ii)i.i'>i be the normalized Gram matrix, or matrix of overlaps, of 
this sample. Consider a subsequence (iVfc) along which the limit in (I2.12p is achieved (now 
with non-random parameters x^) and the distribution of under EG]^^°° converges in 
the sense of convergence of finite dimensional distributions to the distribution of some array 
R°°. For simplicity of notation, let us assume that the sequence {N^) coincides with natural 
numbers. Under KGjf'^°° , the array R^ is weakly exchangeable which means that 

{Rw{l),n{n^ ^ (-^M') {2.18) 

for any permutation vr of finitely many indices. Obviously, this property will be preserved in 
the limit so that R°° is a weakly exchangeable symmetric non-negative definite array and, 
following [6], we will call any such array a Gram-de Finetti array. The Dovbysh-Sudakov 
representation [6] then guarantees that all such arrays are generated by i.i.d. samples from 
random measures on a separable Hilbert space (see [I3j for a detailed proof). 

Proposition 3. If {Ri,v)i,v>i is a Gram-de Finetti array such that Ri i' = 1 then there exists 
a random measure G on the unit ball of a separable Hilbert space such that 

(%')m'>i=(/^'-'^'' + V(1-|I/^'II'))z,.>i (2-19) 
where (p') is an i.i.d. sample from G. 

The importance of the Dovbysh-Sudakov representation in spin glass models was first clearly 
demonstrated in [2] and other examples where this representation played an important role 
can be found in [3], [12] and [16]. Let G be a random measure generating the array R°°, let 
(p') be an i.i.d. sample from G and let Ri^ir = ■ for / ^ V and Ri^i = 1. For simplicity 
of notation, we will now omit oo in R°°. If we denote by (•) the average with respect to G 
then, by (12.171) . the measure G satisfies the Ghirlanda-Guerra identities, 

1 1 " 

ft It 

1=2 

for all p > l,n > 2 and all monomials / of {Ri,i')i^i'<n- Approximating bounded functions of 
the overlaps in by polynomials we also have 

1 1 " 

E(/V^(i?i,„+i)) = -E(/)E(V^(i?i,2)) + - VE(/^/.(i?i,0) (2.21) 
n n ^-^ 

1=2 

for bounded measurable functions / and if). Below, the identities (I2.2ip will allow us to 
identify these asymptotic Gibbs measures but, first, let us show how the limit in (I2.12p can 
be represented in terms of G. By Theorem 2 in [12] , (I2.2ip implies that if q* is the largest point 
in the support of the distribution of Ri^2 under EG"^^ then G is concentrated on the sphere of 
radius i/g* with probability one and, therefore, R is generated by (p' ■ p' + (1 — (f))i,i'>i- 
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Let z{p) and y{p) be two Gaussian processes on the unit ball of our Hilbert space with 
covariances 

mpMp') = ap' ■ p'), mp')yip') = o{p' ■ p'), (2.22) 

let 7] he a standard Gaussian random variable independent of everything else and let 
denote the expectation in r] only. Then the following holds. 

Lemma 2. We have 

Ciqiy/')) (2.23) 

-9iq*)Y/')). (2.24) 

The proof of Lemma[2]is based on the following observation. For a moment, let R = {Ri,i')i,i'>i 
be an arbitrary Gram-de Finetti array such that Ri^i' = 1, let £ be its distribution and let 
G be any random measure generating R as in Proposition [3l It is known that in some sense 
this measure is unique but we will not need it here. Let us define 

$(£) = ElogE,(ch(z(p) + r^(e'(l) - e'dlpin)'/')). (2.25) 

The Gaussian process z{p) here is the same as in (12.221) but we do not assume now that G is 
concentrated on the sphere ||p|p = q*- We will prove that the right hand side in (I2.25P does 
not depend on the choice of the measure G and, indeed, depends only on the distribution C 
in a continuous fashion. 

Lemma 3. The function C — )■ ^{C) defined in Ii2.25\) is continuous with respect to weak 
convergence of the distribution C . 

Proof of Lemma [2l Since R^ is the Gram matrix of the sequence (A^~^/^cr'), we can simply 
think of the measure G~j^ as defined on N~^^'^T,n which is a subset of the sphere ||cr|| = 1 
in M^. Then (12. lip agrees with (12.220 and Lemma [3] implies (12.230 since R^ converges in 
distribution to R°^ and, as we mentioned above, the Ghirlanda-Guerra identities (I2.20p imply 
that G is concentrated on the sphere ||p|P = q* ■ Equation (I2.24p can be proved similarly. 

□ 

Proof of Lemma [3l The proof is almost identical to the proof of Lemma 11 in [15]. For 
simplicity of notation, let us denote 

z,{p) = z{p) + ^{a^) - a\\pr)f" 

and let E^ be the expectation in the randomness of {z{p)) conditionally on all other random 
variables. By standard concentration inequalities for Gaussian processes (see, for example. 
Lemma 3 in [11]) we have that for a > 1, 

P4|logE^^(chz^(p)) -E^logE^(ch^^(p))| > a) < exp(-ca2) (2.26) 

for some small enough constant c that depends only on the function ^ through (12.22p . Since 

< E,logE^(chz^(p)) < log(E,E^chz^(p)) = e'(l)/2, 



lim Elog(chz(<T))_ = E logE^(ch(z(p) + r7(^'(l) - 

N—>oo 



and 



lim Elog(exp?/(cr))_ = ElogE^(exp(|/(p) +77(^^(1) 
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the inequality fl2.26p implies that P(| logE^(chz^(p))| > a) < exp(— ca^) for small c and large 
enough a and, therefore, if we denote log^o; = max(— a, min(logx, a)) then for large a, 

|ElogE^(ch;z^(p)) -Elog„E^^(ch2^(p))| < exp{-ca^). (2.27) 

Next, if we define chaX = min(chx,cha) then using that | log^x — log^y\ < e°'\x — y\ and 
|chx — ch(ix| < chx/(|a;| > a) we can write 

|Elog„E^(chz^(p)) - Elog„E^(ch„z^(p))| < e''E(|ch2;^(p) - chaZ^(p)|) 

< e'^E(chz^(p) J(|z^(p)| >a)>. 

By Holder's inequality, this can be bounded by 

e'^(E(E,,,ch2z^(p)))'/' (E(P,,,(|z,(p)| > a)))'^' < exp(-ca2) 

for small c and large enough a since P^_^(|z^(p)| > a) < exp(— ca^). Combining with (12.271) . 

|ElogE^(ch2^(p)) -Elog„E^(cha2^(p))| < exp(-ca2). (2.28) 

Approximating logarithm by polynomials on the interval [e~'^,e"], E log^ E^(cha2;^(p)) can 
be approximated by a linear combinations of moments 

E(E,(ch,2;^(p)))" =E(E,E, J]ch,(2;,KpO)> (2-29) 

l<r 

where we used replicas and where z^i{p^) = z{p^) +rf{^'{l) — ^'(Hp'IP))^^^ and (?]') are i.i.d. 
standard Gaussian. Since the covariance of the Gaussian sequence [zy^i[p^)) is equal to 

r(p' • /) + 5M'(r(l) - e'dlp'f )) = ^M'), 
the function inside the Gibbs average on the right hand side of (I2.29P is equal to 

E,E^ W dia{z^i (p')) = F{{R^,yi,<r) (2.30) 

l<r 

for some continuous bounded function F of the overlaps {Ri,i')i,i'<r- Together with (12.281) this 
shows that we can approximate ^{C) arbitrarily well by a linear combination of K{F{R)) 
for some continuous bounded functions F of finitely many overlaps, which proves that ^{C) 
is continuous with respect to the distribution £ of the overlap array R. 

□ 

Identifying asymptotic Gibbs' measures using ultrametricity. To show that the 
lower bound in (I2.12p matches Guerra's upper bound, it remains to identify the difference of 
(I2.23P and (I2.24p with the functional (II. 7p . Since the asymptotic Gibbs measure G satisfies 
the Ghirlanda-Guerra identities (I2.2ip . the main result in [T7] implies that the support of G 
is ultrametric with probability one, that is, 

E(j(i?i,2 > min(i?i,3, i?2,3)) > = 1. (2.31) 
9 



Given k > 1, let us consider a function ^(g) on [0, 1] such that 

K{q) = j/k for j/k <q<{j + l)/k, j = 0,...,k. (2.32) 

Equation fl2.3ip imphes that for any q the inequahty q < ■ defines an equivalence 
relation / ~ /' and, therefore, the array (/(g < Ri^i'))i,i'>i is non-negative definite since 
it is block-diagonal with blocks consisting of all elements equal to one. This implies that 
= {k,{Ri,i>))i.i>>i is non-negative definite since it can be written as a convex combination 

In addition, it is clear that is weakly exchangeable and satisfies the Ghirlanda-Guerra 
identities fl2.2ip . Then, by the Dovbysh-Sudakov representation (12.191) . can be generated 
by a sample from some random measure G"* on the unit ball of a Hilbert space. If for 
simplicity we assume that q* ^ j/k for j < k then K{q*) is the largest point in the support 
of the distribution of k{Ri^2) and by Theorem 2 in [12] the measure is concentrated on 
the sphere ||p||^ = ^(g*). When A; — )• oo, the distribution of R'^ converges weakly to the 
distribution of R and if we denote by the average with respect to the measure G"* then 
Lemma [3] implies that 

ElogE,{ch{z{p)+v{C'{l)-C'W)y^'))^ (2.33) 

approximates the right hand side of (12.231) . It is well known that an ultrametric measure, 
such as C^, that satisfies the Ghirlanda-Guerra identities and under which the overlaps R'^ 
take finitely many values as in (I2.32p can be identified with the discrete Ruelle probability 
cascades by the Baffioni-Rosati theorem [4] (see the proof of Theorem 15.3.6 in [25] for 
details). The fact that in this case (I2.33P coincides with Xo{m,q) in (11.70 with parameters 
qj = j/k and rrij = E(/(i?i^2 < Qj)) is also well known (see e.g. Theorem 14.2.1 in t25i|). One 
can similarly show that (12.240 corresponds to the second term in (11.70 and this finishes the 
proof of Theorem [H One could also work with the continuous Ruelle probability cascades 
using a general theory developed in [5] but, at this point, it was easier to simply discretize 
the overlap array and Gibbs measure. 

□ 
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